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In quantum systems, one usually seeks to minimize dephasing noise and disorder. The efficiency of 
transport in a quantum system is usually degraded by the presence of noise and disorder. However, 
it has been shown that the combination of the two can lead to significantly more efficiency than 
each by itself. Here, we consider how the addition of nonlocal noise, in the form of incoherent 
hopping, affects the transport efficiency. We show that incoherent hopping introduces additional 
local extrema in the efficiency function and investigate the crossover from a quantum random walk 
to a classical random walk. 
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I. INTRODUCTION 

Environmental noise is simultaneously intrinsic to any 
physical process and usually detrimental to its efficient 
operation. In the context of excitation transport in 
quantum systems, dephasing noise leads to decoherence 
and reduces the transport efficiency. However, in sys¬ 
tems possessing spatial or dynamical symmetry break¬ 
ing, fluctuating noise can at times enhance transport, 
as in ratchets and Brownian motors [I|. The idea of 
stochastic resonance has been used to describe transport 
and power enhancement in a variety of classical, biolog¬ 
ical and quantum networks 0. Among those, a natural 
light-harvesting system exhibiting extremely high energy 
transfer efficiency has attracted a lot of attentions and 
even helped precipitate the field of quantum biology Q. 
The mechanism of excitation energy transfer (EET) in 
light harvesting complexes and to what extent EET is 
influenced by quantum mechanical effects, calls for a con¬ 
sistent theoretical model Q. 

Recently, it has been pointed out that for realistic cou¬ 
pling strength in photosynthetic pigments, quantum and 
classical coherent transports are identical . Previously, 
a quantum walk @ framework was already presented 
for disordered systems, in which under the influence of 
quantum interference, an interplay between free Hamil¬ 
tonian evolution and thermal fluctuation in the environ¬ 
ment leads to an increase in EET. In the absence of noise, 
there is hardly any transport since disorder leads to An¬ 
derson localization of eigenstates. However, the addition 
of a local dephasing noise scrambles the wave function so 
that it is no longer localized in space, leading to increased 
transport efficiency Q ■ In Ref. 0 , the dephasing-assisted 
transport in quantum networks and biomolecules was in¬ 
vestigated assuming a similar system-environment inter¬ 
action. They showed that, reminiscent to stochastic res¬ 
onance in classical system Q, transport of excitations 
across dissipative quantum networks can be enhanced by 
local dephasing noise. 

In this work, we study the effect of nonlocal noise, in 
the form of incoherent hopping, on the transport effi¬ 


ciency in a prototypical one dimensional spin chain. We 
obtain analytical and numerical results for transport ef¬ 
ficiency in quantum chains of varying sizes. We show 
that in the absence of disorder, the addition of incoher¬ 
ent hopping leads to a local minimum in the efficiency, 
marking the crossover from coherent quantum to inco¬ 
herent classical dynamics. We show that in a disordered 
chain, incoherent hopping leads to a local minimum and 
a local maximum. We also explore an experimental im¬ 
plementation of such noise-enhanced transport dynamics 
with trapped ions. 


II. GENERAL FORMALISM 


We adopt the quantum walk framework of Refs. [1, Q 
to a one-dimensional chain of N sites subject to dissipa¬ 
tion. Each site is a spin-1/2 with an energy difference 
between the two spin states. The Hamiltonian and the 
Lindblad superoperators are 
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where {kl) denotes summing over nearest neighbors and 
{A, B} = AB-\-BA. There is coherent hopping v between 
nearest neighbors. Each site dissipates energy into the 
environment at a rate T when it is in the upper energy 
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state, corresponding to the decay of an excitation. In 
addition to the local dephasing rate 7 already considered 
in previous studies [^, , we also include an incoherent 

hopping rate y/j. This rate is so-called since it originates 
from the fluctuations in the nearest-neighbor hopping. 
We assume white-noise properties in these exterior field 
so that the system undergoes a Markov evolution as in 
Ref. [13 . 

To study the energy transfer, we add an additional site 
TV -I-1, which is dissipatively coupled to site N with a rate 
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Site TV -I- 1 is a sink as once the excitation enters, it is 
removed from the system. 

The Lindblad master equation for the density matrix 
p reads 

^ [H, p] + Cdissip) + Cdephip) + £sink(p), (5) 

This is basically a set of TV^ linear ordinary differential 
equations with initial conditions po and thus can be con¬ 
veniently solved with the Laplace transformation. 

We assume that the chain is initially excited at site 1. 
The total population depends only on T and A, while the 
relative population between sites depends on the other 
parameters. Throughout this work, we choose A/T = 
10.0 and T/v = 0 . 02 . 

Transfer efficiency is defined as the probability of ex¬ 
citations arriving at the sink, 

poo 

r] = A p]sf,Nir)dT. ( 6 ) 
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Below, we explore how the inclusion of incoherent hop¬ 
ping affects the efficiency of transfer. 


III. UNIFORM SYSTEM 

We start with a chain of TV = 2 sites. In a rotating 
frame, the coherent dynamics of this system is governed 
by two parameters, the coherent hopping v and the dif¬ 
ference in site energies 6 = uJi — U! 2 - It is well-known 
that in a closed quantum system, coherent hopping fa¬ 
vors excitation transfer while disorder tends to localize 
excitations [HI- Therefore, the closed quantum system 
has two inherent competing energy transfer processes. It 
has been shown that one can counteract the localizing 
effect of disorder by including dephasing 7 . 

Using Laplace transformation, the EET efficiency can 


be obtained analytically, 

r(r -h A) 

r(r + A) + ( 4 r; 2 /A + 7?.)(2r + A)J ’ 
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A = v+—, (8) 

P = 2r + A + 2(7 + 7 ^). (9) 

Here, V is the sum of all incoherent rates. 

First, we set disorder (5 = 0 to see the competition 
between coherent hopping v and incoherent hopping jh- 
Since local dephasing 7 only shifts the total dissipation, 
for simplicity, we also set it to be zero. Straightforward 
algebra shows that for a given set of T and A, the effi¬ 
ciency as a function of 'jh/v has a local minimum, as seen 
in Fig. [TJa). The reason for the minimum is as follows. 
When 7 /i is zero, the dynamics is fully coherent, and 
the excitation undergoes a quantum random walk. As 
7 /j increases, the efficiency decreases since the incoherent 
hopping destroys the coherence required for a quantum 
random walk. However, when ■jh is sufficiently large, the 
dynamics is dominated by incoherent hopping and the 
excitation executes a classical random walk. The effi¬ 
ciency minimum marks the transition from a quantum 
to a classical random walk. 

To characterize the local minimum, we define the depth 
of the dip as the difference between the minimum 77niin 
and the saturated efficiency rjc when either coherent or 
incoherent hopping is much stronger than dissipation. 


A 

2r-k A 



FIG. 1. (Color online) Competition of coherent and incoher¬ 
ent hopping for a uniform (<5 = 0) two-site chain, (a) Energy 
transfer efficiency as a function of ratio ■yh/v for A/T = 10 and 
r/a = 0.02; The minimum near yn = v marks the transition 
from a quantum to a classical random walk, (b) Dependence 
of dip position on total dissipation rate; (c) Dependence of 
relative dip depth on A/F for T/v = 0.02; and (d) Depen- 
dend of relative dip depth on (2r -|- A)/v for A/F = 10. Red 
solid lines are for zero local dephasing and blue broken for 
y/v = 0 . 2 . 
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FIG. 2. Energy transfer efficiency as a function of ^hjv for a 
uniform (<5 = 0) chain of various length N. Local dephasing 
7 is taken to be zero. For all the plots, A/F = 10 and V jv = 
0 . 02 . 


Consequently, the relative depth is 
At? ^ 7?e - ^7min 
?7c ?7c 

=_ piAA_ ^ ( 10 ) 

r(r + A) + [2V2v - (2r + A)] (2r + A) 

Figures [Hb-d) illustrate the dependence of dip posi¬ 
tion and relative depth on two different ratios A/F and 
(2r + A)/t;. In the same figures, we include the transfer 
efficiency with a nonzero local dephasing ^ jv = 0.2 (blue 
lines) for comparison. Incorporating local dephasing 7 
shifts the minimum to the left and could make the de¬ 
creasing part of the curve invisible. At the same time, it 
also makes the contrast between the dip and saturation 
more evident. 

For longer chains, analytic results similar to Eq. 0 
can in principle be obtained. However, the tedious closed 
form results are not illustrative. Instead, we resort to 
the numerical evolution of the master equation. Figure [2] 
presents the simulation for chains of iV = 3, 5, 10, and 20 
sites respectively. As N increases, the efficiency decreases 
while the relative depth of the minimum increases. 


IV. DISORDERED SYSTEM 

Now we turn to systems whereby disorder (5 plays a 
role. It is seen from Eq. o that 7? is a monotonic in¬ 
creasing function of the combined parameter 

f = ( 11 ) 

Simple algebra shows that depending on the ratio S/v, 
the efficiency as a function of ‘jh/v is basically catego¬ 
rized to two different types: (1) when 6/v > 1/2 [region 



FIG. 3. (Golor online) (a) and (b) Parameter space partition 
for different types of transfer efficiency as a function of the 
ratio of incoherent vs coherent hopping for a generic two- 
site {N = 2) chain. In the lower right region oi v — S space 
denoted by B, efficiency monotonically increases with 'yh/v as 
seen in (d). In contrast, in the upper left region A, for suitable 
combination of other incoherent rates [in (b)], efficiency may 
show local extrema by varying 7 ti/u. In (b), dissipation rate 
2r +A and disorder 5 space is divided by T>± into region I, II, 
and III where efficiency displays two, one, and no extrema, 
(c) and (d) Efficiency as a function of yh/v for two typical 
parameter sets. In both subplots, A/F = 10.0, T/v = 0.02, 
red (solid) lines for 7 = 0 and blue (broken) lines for y/v = 
0 . 2 . 


B in Fig. [HKa)], the efficiency curve does not have local 
minimum or maximum; (2) when S/v < 1/2 [region A 
in Fig. EKa)], there may be a local maximum potentially 
lurking in the curve in addition to the aforementioned 
local minimum, although the existence of local extrema 
will depend on other parameters. 

Figure [3| shows the parameter space partition for dif¬ 
ferent types of transfer efficiency as a function of yn/v 
and the characterization of the extrema for a two-site 
{N = 2) chain. In Fig. [31(a), the lower right region of 
the v-S space is the territory where efficiency displays no 
local maximum. Figure [3Kb) further delineates three re¬ 
gions I, II, and III when disorder is not large compared 
to the coherent hopping {6/v < 0.5), where efficiency as 
a function of yn/v shows both local extrema, only local 
minimum, and no local extrema, respectively. As before, 
the blue broken line illustrates that the nonzero local de¬ 
phasing shrinks the region where extrema actually show 
up. In Fig. He) we deliberately choose S/v = 0.4 and 
7 = 0, which falls into the region I of Fig. Hb), and plot 
the efficiency as a function of yt/v. As expected, both 
extrema show up. FigureHd) presents another efficiency 
curve for the parameter set in the region B of Fig. Ha). 

The incoherent hopping yh has been neglected in pre¬ 
vious studies, on the assumption that it is much smaller 
than the local dephasing 7. We find that a small ra¬ 
tio yh/'y does lead to a visible competition between the 
coherent and incoherent hopping mechanisms. Figured) 
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FIG. 4. (Color online) Energy transfer efficiency as a function 
of local dephasing rate as from analytical formula Eq. 0 for 
various ratio ^ of hopping fluctuation to local dephasing rate 
in disordered S/v = 5.0 two-site system, where A/E = 10.0 
and T/v = 0.02. 


FIG. 5. Energy transfer efficiency in a three-site Lambda-type 
chain as a function of (a) ratio of disorder to cohorent hopping 
for = 7 = 0; (b) ratio of incoherent hopping to coherent 
hopping for 7 = 0; (c) ratio of local dephasing to coherent 
hopping for yh = O.OOI7; and (d) ratio of local dephasing 
to coherent hopping for jh = 0. For (b)-(d), A/F = 10, 
r/v = 0.02, and 5/v = 10. 


plots the EET efficiency as a function of y/u for various 
'jhll- The dips are the turning points where “dressed” 
coherent dynamics totally gives way to incoherent hop¬ 
ping dynamics, such that the system can effectively be 
described by a classical master equation instead of the 
Lindblad master equation. 

From Eq. (HU, we infer that when the two additive 
terms are equal, the system is balanced by two separate 
influences, dressed coherent dynamics and purely inco¬ 
herent dynamics. Therefore by solving = 7/1, we 

obtain the critical incoherent hopping 7^. We find that 
there exists a practical parameter region where 7^ reaches 
the smallest scale as common belief states. The combined 
parameter / can thus be seen as an effective hopping rate 
due to the competition among various coherent and in¬ 
coherent processes within the one-excitation manifold. 

Figure [5] displays results for a disordered Lambda-type 
three-site chain. We see that a nonzero saturated trans¬ 
fer efficiency rate rjc, while in the absence of incoherent 
hopping, strong enough local dephasing leads to zero ef¬ 
ficiency. 

V. EXPERIMENTAL REALIZATION 

Throughout this work, we have scaled system disor¬ 
der and all incoherent rates with respect to the coher¬ 
ent interaction between neighboring spins. To make con¬ 
tact with experimentally realizable conditions, we need 
to know the typical coherent interaction strength. In 
a recent work by Z. Meir et al. [T^, coherent coupling 
between distant Sr+ ions was attributed to the far-held 
resonant dipole-dipole interaction. They performed a di¬ 
rect spectroscopic measurement of the cooperative Lamb 
shift at the 55'i/2 -O- 5Pi/2 optical dipole transition fre¬ 


quency. The shift can be expressed as the expectation 
value of the interaction in the excited state. For two 
ions, the cooperative shift was predicted to be 130 kHz 
at a distance of about 5 fxm. However, in other ion-trap 
experiments, depending on the specihc ion species and 
different trapping techniques, the coherent dipole-dipole 
interaction can be dramatically different. Recently, two 
ion-rap experiments aimed to produce tunable long-range 
coupling with ^^^Yb+ 0 and ^°Ca+ 0. The for¬ 
mer experiment utilized the ‘^Si/ 2 \F = 0, = 0) and 

'^Fi/ 2 \F = = 0) hyperhne ‘clock’ states in a mag¬ 

netic dipole transition, while the latter experiment used 
the \Si/ 2 ,'m = +1/2) and |Zl5/2,m- = +5/2) in an elec¬ 
tric quadrupole transition. By adjusting a combination 
of trap and laser parameters, both groups can vary the in¬ 
teraction range in addition to actually providing nearest- 
neighbor interaction, on the order of 0.1 to 1 kHz. In 
Ref. [3 1 it is pointed out that electric field noise leads to 
imperfect conservation of excitation number, while laser- 
frequency and magnetic-field fluctuations give rise to de¬ 
phasing during the dynamics in the linear Paul trap. 


VI. CONCLUSION 

We have studied excitation transfer in a one¬ 
dimensional chain subject to both local dephasing noise 
and incoherent hoppping. An analytical form of the en¬ 
ergy transfer efficiency is given for a two-site chain, which 
predicts rich behavior due to the competition between 
different processes. The main features also appear in 
longer chain. The phenomena predicted can be tested in 
quantum many-body system of trapped atomic ions re- 
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cently reported to exhibit cooperative Lamb shift [T^ or 
tunable for long-range interactions [H, [3 • 

The quantum walks model we have based this work on 
adopts the Born-Markov approximation and deals with 
the dynamics to second order in the sj^em-bath cou¬ 
pling. There have been several groups |15l-[l^ trying to 
develop fully non-perturbative and non-Markovian ap¬ 
proaches for the open system quantum dynamics. We 


expect the non-local fluctuation introduced in this work 
to be an indispensable part of the more generic theory. 
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